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Abstract 



In this paper we investigate the possible existence of multi-point correlation in the 
arrival direction of the UHECR events detected by the High Resolution Fly's Eye 
(HiRes) stereo detector. Multi-point correlations could result from the deflection 
of UHECRs by galactic and intergalactic magnetic fields, and the subsequent dis- 
persion of arrival directions from point like sources. The search is performed by 
calculating the solid angle subtended by the polygon between triplets and quadru- 
plets of events in the HiRes data. The resulting distribution of solid angles is then 
compared to the cumulative distributions from multiple simulated isotropic data 
sets to estimate the significance of any excess. We also looked for potential correla- 
tion of the small solid angle triplets found in the data with the locations of the BL 
Lacertae (BL Lac) objects. Neither statistically significant clustering nor significant 
correlations with BL Lac objects were found in these studies. 

Key words: Anisotropy, Ultra High Energy Cosmic Rays, HiRes, Flys Eye, Air 
Fluorescence, BL Lacertae. 



1 Introduction 



While there have been reports of possible correlation, the sources of UHECRs 
are currently unknown ([1], [2], [3]). One of the challenges in answering this 
question lies within our limited knowledge of the galactic, inter- and extra- 
galactic magnetic fields. Deflections of UHECR particles in the magnetic fields 
complicates the correlation that should otherwise exist between the arrival di- 
rection of these particles and the source position in the sky. However, moderate 
deflections of particles from the same source through the magnetic field should 
result in clustering in their arrival direction. 

Clustering in the arrival direction of UHECRs using the HiRes stereo data 
set has been investigated in a previous study searching for events in the sky 
with an arrival direction lying on a great circle [4]. In this paper, however, the 
clustering search uses the solid angle subtended by polygons in the data. This 
method does not look for a particular structure (e.g. an arc). The search using 
this method could involve a coherent or an incoherent clustering of events. 

In addition to looking for events clustering, a search for correlation with BL 
Lac objects is performed. BL Lacs are blazars, a type of Active Glactic Nu- 
clei (AGN) with its jets aligned along our line of sight. Blazars are an estab- 
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lished source of TeV Gamma rays [5] and they are considered highly motivated 
sources for UHECRs. 

Recent studies ([6], [7], [8], [9]) have reported correlations between the AGASA, 
Yakutsk, and HiRes cosmic ray data sets and a subset of BL Lac objects. In 
addition, a study in [1] verified the correlation between the HiRes stereo data 
set and a subset of BL Lac objects (BL Lacs with magnitudes < 18, includ- 
ing the Highly Polarized [HP] BL Lacs). However, and as mentioned in [1] 
the significance of the observation is difficult to estimate due to the apparent 
a posteriori nature of the initial search. These reports have yet to be confirmed 
with a clean, independent data set. 

In this paper we apply alternative methods to examine this HiRes stereo data 
set for three- and four-point correlations as well as possible correlation between 
clusters found in the data and BL Lac objects. The data set used in this paper 
is the High Resolution Flys Eye (HiRes) stereo data set collected between 
12/1999 and 1/2004 and consists of 271 events with energy > 10 19 eV. This 
is the same data set used in ([1], [10]). Particles with lower energies would 
likely suffer deflections large enough to destroy the correlation between the 
arrival direction and source position. The angular resolution of the HiRes 
stereo detector for this data set is ~ 0.6°. More details about this data are 
given in [10]. 



2 Search for Three-point correlation 

In order to identify statistically significant clustering, we look for excess in 
the real data over what is found in simulated isotropic background sets. The 
simulation data sets used for this analysis are based on a Monte Carlo that ac- 
counts for full geometrical and exposure acceptance of the detector. A detailed 
description of this MC can be found in [10]. 

We examined every possible combination of three events in the data set. How- 
ever, because we are looking for clustering around a point source, we con- 
centrate on those triplets subtending small solid angles. Small solid angles 
formation could imply the presence of an arc formation searched for in [4] 
and/or clustering in the arrival direction of detected UHECRs. To this end, 
we restrict our search to have angular separation between any pair of events 
within a triplet of less than 10°. The solid angle f2, is calculated using Girard's 
theorem ([11], [12]). 

Figure 1 shows the frequency with which groupings of events appear with a 
given spherical solid angle. The data are shown as (+) points while the solid 
lines mark the mean and the ± la limits from 1000 simulated data sets. The 
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data lies well within the isotropic MC expectations, with no excess is seen at 
small solid angles. Note that the number of points (the data points) that lie 
outside the ± la are significantly less than 32% of the total number of points. 
This is due to the bin-to-bin correlation which we have studied and reported 
in [13]. The chance probability for the data to have at least as many triplets 
in the smallest angle bin (smallest Q) is calculated by counting the number 
of times the simulated data sets yield a stronger signal than the data. The 
chance probability of the signal to be greater than 14 (the value of the first 
bin from the HiRes data) is found to be equal to 0.75. Hence, no significant 
signal is found for small scale clustering of the HiRes stereo data set using the 
three-point correlation method. 




Q (steradians) 



Fig. 1. The distribution of solid angles, f2, of the spherical triangles. The points (+) 
represent the spherical triangles found in HiRes stereo data set. The lines represent 
the mean and the ± la limits for the simulated data sets. The opening angles 
between each pair of points within the three points are required to be < 10°. 

3 Three Point Correlation With BL Lac Objects 

As an alternate approach to find a possible correlation between the HiRes 
stereo data with the BL Lac objects, we used a modified three-point correla- 
tion method. A primary motivation for performing a multi-point correlation 
search with BL Lacs was that the two-point correlation search reported in [1] 
was consistent with the experimental angular resolution of the HiRes detector 
without any magnetic smearing. A possible interpretation of this result was 
that the events correlated to the BL Lac objects are from neutral primaries, 
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which is inconsistent with HiRes's report of proton-dominance of UHECR 
composition at the highest energies [14]. In the three- and four-point cor- 
relation studies we test the more pedestrian hypothesis that the two-point 
correlation seen might simply be part of a more smeared out multi-event cor- 
relation. 

In this study we calculated the solid angle of the spherical triangles contained 
between two event points from the data and one point from the BL Lac set. 
As in the case of the self three-point correlation study presented above, we 
restricted our search to those events with angular separation less than 10°. 
Figure 2 shows the distribution of solid angles from the spherical triangles 
found in the real data set shown as points (+) to be compared to the mean and 
± la of 1000 simulated data sets shown as solid lines. As would be expected 
from the observation of a two-point correlation with BL Lacs. In the plot, we 
do see an apparent excess at small solid angles. Here the chance probability 
for the value of the real data to be greater than 55 (the value of the first bin) 
is found to be equal to 5.4% (~ 2a). 

It is important to note that the significance of the three-point correlation is 
much weaker than that found in [1] for the two-point correlation, which is 
consistent with what is expected from a random third point association with 
the existing correlated pairs. Therefore, no new information was found on BL 
Lac correlation with this study. 



4 Cross Correlation Between the Centroid of the Triplets With 
Small Solid Angles and BL Lac Objects 

In addition to the test in the previous section. We searched for cross-correlation 
between the small solid angle clusters found in our data and the BL Lac 
objects. The cross-correlation is tested for triplets of solid angle Q < 0.004 
sr. This cut corresponds to the selection of triplets that lie in an arc of up to 
10° and allows for a deflection of 1° for the case where the middle event is 
centered between the two end events (consistent with the angular resolution of 
HiRes). Figure 3 illustrates such a triplet. A total of 207 triplets were selected 
for study by this criterion. The comparison was performed for the same set of 
BL Lacs used above. 
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Fig. 2. The distribution of areas for the spherical triangles found in using BL Lacs 
as the third point. The points (+) represent the spherical triangles using the HiRes 
stereo data set with the BL Lacs, where two points are from the data set and the 
third point is from the BL Lac location. The lines represent the mean and the ± la 
for the simulated data sets with the same study. The opening angles between each 
two points of the three points are required to be < 10°. 




10° 

Fig. 3. A diagram of a three points with a separation angle between the furthest 
two points of 10°. The middle point of the triangle is displaced one degree up and 
the end two points are displaced one degree (consistent with angular resolution 
smearing) down from where the three points would be placed on the same arc. 
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Figure 4 shows the histogram of cos(8), where 9 is the opening angle between 
the centroid point of a spherical triangles with solid angle, Q < 0.004 sr and 
the BL Lac set. The points (+) represent the data set and the solid lines show 
the mean and ±lcr bounds from 1000 isotropically generated MC sets. As a 
figure of merit, we used the excess significance of the last bin where cos(9) is 
closest to 1 (9 of < 10°). Looking at the 1000 isotropic simulated data sets, we 
found the chance probability to be 0.3 of having an excess at least this great. 

Note that a few bins in this histogram display an excess of approximately 2a 
in their occupancy level. However, one must remember that the values of the 
bins in this distribution are correlated. Therefore, each bin can not be treated 
independently. Consequently, the values of these bins can not be combined. 
Any excess seen is therefore no more than 2a and consistent with the three- 
point correlation in the previous section. For this study we made the a priori 
choice to use the significance of the last bin and that is what we consistently 
use. Consequently, we conclude that there is no significant excess between the 
small solid angle triplets with the BL Lac objects. 




Fig. 4. The frequency distribution of cos(9), where 8 is the opening angle between 
the middle point of the spherical triangles and the BL Lac. The spherical triangles 
have area, f2 < 0.004 sr and the BL Lacs have magnitude < 18. The points (+) 
represent the real data set. The lines represent the mean and ±1<7 limits calculated 
from 1000 isotropic MC sets. 
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5 Four Point Correlation 



In the previous correlation study we looked for three-point correlation in the 
HiRes stereo data set using their spherical triangle area. In this section we 
will attempt to extend the technique to four-point correlations in an attempt 
to further investigate potential clustering in the HiRes stereo data set. We do 
this by calculating the solid-angle bounded by each quadruple set of events in 
the data. The enclosed solid angle is calculated by dividing the quadrilateral 
into two spherical triangles. The actual division can be done in two different 
ways, as illustrated in Figure 5. The total solid angle is simply the sum of 
these two spherical triangles. Note that the value of the total solid angle is the 
same for either of the two divisions. In addition, the events can be configured 
such that a point of the quadruple is inside a spherical triangle formed by 
the other three points as shown in Figure 5. In such the solid angle is 

that of the spherical triangle bounded by the outer three points. The resulting 
histogram of the solid angle areas is then compared to simulated data sets to 
estimate the significance of any excess. 



Fig. 5. The first two plots describe two different divisions of a spherical quadrilateral 
into two spherical triangles. The same overall solid- angle is obtained from the sum 
of the two triangles, regardless of which division is used. The last plot shows the case 
where one of the points of quadruplet lies inside a spherical triangle bounded by 
the other three. In such case the solid-angle is obtained from the spherical triangle 
containing the fourth point. 

As in the case of the three-point autocorrelation study, we compare the dis- 
tribution of the solid angles from the real data to that found in the simulated 
isotropic MC sets. The only difference is that we are using quadruplets rather 
than triplets. The solid angle of the quadruplets is calculated as the area of 
the exterior spherical triangles [13]. As in the case of the three point study, 
we have made a cut on the angular separation between pairs of points in the 
quadruples to be less than 10°. 

Figure 6 shows the occupancy level vs. the solid angle area of quadruples 
with the 10° angular separation cut. As previously, the actual HiRes data 
set is shown by the data points. The mean and ± ler intervals are shown 
by the solid lines. The figure clearly shows no obvious excess in the lowest 
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bins. Again, the chance probability of the first bin from 1000 MC sets area 
(corresponding to the smallest solid angle) is again calculated by dividing the 
number of simulated data sets that yielded a stronger signal than that of the 
real data set by the total number of simulated data set. In this case, the value 
was found to be ~ 0.70. Hence, clearly no significant four-point clustering is 
found in the data. 



When looking at the distribution of data points in Figure 6, they are found 
to be systematically lower when compared with the mean obtained from the 
simulated sets. This deficit suggests that our understanding of the detector 
aperture is not good enough to reproduce the data beyond a four-point cor- 
relation search. In addition, the number of clusters of small angle four-point 
cluster is barely significant for statistical analysis. Hence, we decided in this 
study to stop at four-point correlation and did not proceed to five (or more) 
point auto correlation. 
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Fig. 6. The histogram of the solid angle between each of four points with the angular 
separation cut. The points (+) represents the real data sets, while the solid lines 
represent the mean and the ± la of a 1000 simulated data set. 



6 Conclusion 



In this paper we described a search for small scale clustering in the HiRes stereo 
data set using three point correlation. The search is applied by histogramming 
the area subtended by spherical triangle defined by triplets in the HiRes stereo 
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data set. We required that every pair in a triplet lie within 10° of angular 
separation. The result from the histogram is then compared to that found 
in the simulated data sets. The simulated data sets were generated to have 
the same exposure as that of the real observed data. No excess of small scale 
clustering in the HiRes stereo data set was found. 

We also looked for correlation between the HiRes stereo data set and the BL 
Lac objects using the three point correlation technique. The strategy used in 
this study is looking for three points where two of the three points are from the 
data set and the third one is a BL Lac location, each of which has a maximum 
angular separation of 10°. The significance found is approximately 2a and is 
weaker then the the two-point correlation study in [1]. We conclude that no 
new information on BL Lac correlation was obtained from these studies. In 
addition, we looked for cross correlation between the HiRes stereo data triplets 
with £7 < 0.004 sr and the BL Lac set. No significant correlation was found 
between the triplets and the BL Lac objects. 

Finally, we tested for four point correlation in the HiRes stereo data set and 
found no statistically significance excess. We conclude that the four point auto 
correlation study is pushing both statistical borderline and systematic limits, 
of the reliability of the simulation. For this reason, we decided not to extend 
the study to four point BL Lac correlation or to five point auto correlation. 
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